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We reveal the center vortex content of SU{2) calorons and ensembles of them. We use Laplacian 
Center Gauge as well as Maximal Center Gauges to show that the vortex in a single caloron consists 
of two parts. The first one connects the constituent dyons of the caloron (which are monopoles in 
Laplacian Abclian Gauge) and extends in time. The second part is predominantly spatial, encloses 
one of the dyons and can be related to the twist in the caloron gauge field. This part depends strongly 
on the caloron holonomy and degenerates to a plane when the holonomy is maximally nontrivial, 
i.e. when the asymptotic Polyakov loop is traceless. Correspondingly, we find the spatial vortices in 
caloron ensembles to percolate in this case. This finding fits perfectly in the confinement scenario 
of vortices and shows that calorons are suitable to facilitate the vortex confinement mechanism. 



I. INTRODUCTION 

To answer the question of what drives confinement 
and other nonperturbative effects in QCD, basically three 
sorts of topological excitations have been intensively ex- 
amined over the years; instantons, magnetic monopoles 
and center vortices. Instantons as solutions of the equa- 
tions of motion are special: they are the relevant ob- 
jects in a semiclassical approach. While the generation 
of a chiral condensate is very natural via the (quasi) 
zero modes ^, confinement remained unexplained in this 
model. 

At finite temperature, where the classical solutions are 
called calorons [1-3], there has been quite some progress 
recently, due to two effects. First of all, the asymptotic 
Polyakov loop plays a key role determining the proper- 
ties of a new type of caloron solutions (for more details 
on calorons see Sect. II and the reviews [4]). Under the 
conjecture that the asymptotic Polyakov loop is related 
to the average Polyakov loop, the order parameter of con- 
finement, calorons are sensitive to the phase of QCD un- 
der consideration. 

Secondly, the new calorons with nontrivial holonomy 
consist of N dyons/magnetic monopoles for the gauge 
group SU{N). In this way, contact is seemingly made 
to the Dual Superconductor scenario. We stress that the 
dyon constituents of calorons appear in an unambiguous 
way as classical objects. 

This is in contrast to Abelian monopoles and also cen- 
ter vortices, the other sorts of objects used to explain 
confinement. In the sense, that they are widely accepted, 
they are not of semiclassical nature. They represent 
gauge defects of codimension 3 and 2, respectively, which 
remain after the respective gauge fixing and projection. 



^ This mechanism is based on the index theorem and thus will 
work; for any object with topological charge. 



Their interrelation and the fact that they are a prerequi- 
site for the occurrence of topological charge in general 
has been quantitatively studied in the past [5]. 

Monopoles are usually obtained by applying the Max- 
imal Abelian Gauge (MAG). Center vortices in lattice 
QCD can be defined through a center projection (there- 
fore called P-vortices) after the lattice gauge field has 
been transformed into the Maximal Center Gauge (di- 
rectly by Direct Maximal Center Gauge [DMCG] or in- 
directly by Indirect Maximal Center Gauge [IMCG] , with 
the MAG as preconditioner) or into the Laplacian Center 
Gauge (LCG). We refer to Sect. Ill for the technicalities. 

We would also like to mention another important de- 
velopment during recent years. Fermionic methods have 
become available to study topological structures with- 
out the necessity of smoothing. Singular, codimen- 
sion 1 sheets of sign-coherent topological charge have 
been found and proposed to be characteristic for gen- 
uine quantum configurations [6, 7] and potentially im- 
portant for the confinement property. The relation to 
the other low-dimensional singular topological excita- 
tions is still not completely understood. In this scenario, 
(anti)selfdual objects like calorons typically appear as 
topological lumps after smearing [8], i.e. at a resolution 
length bigger than the lattice spacing [7]. 

The physical mechanisms assigned to calorons would 
be based on their quantum weight [9] , their moduli space 
metric [10], the particular features of their fermionic zero 
modes [11] and the specific suppression of dyons accord- 
ing to the action they acquire in different phases [12] . The 
last observation suggests an overall description of confine- 
ment [13] and deconfinement in terms of calorons' dyon 
constituents as independent degrees of freedom. The pro- 
posed generalized (approximative) moduli spaee metric, 
however, presents some difficulties [14] which are not yet 
overcome. 

Since calorons unify instanton and monopoles, it is 

natural to ask for the relation to center vortices. In 
four dimensions the latter are two-dimensional world- 
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sheets with Wilson loops taking values in the center of 
the gauge group if they are linked with the vortex sheet. 
Vortices that randomly penetrate a given Wilson loop 
very naturally give rise to an area law. Since vortices are 
closed surfaces, the necessary randomness can be facili- 
tated only by large vortices. This is further translated 
into the percolation of vortices, meaning that the size 
of the (largest) vortex clusters becomes comparable to 
the extension of the space itself. This percolation has 
been observed in lattice simulations of the confined phase 
[15], while in the deconfined phase the vortices align in 
the timelike direction and the percolation mechanism re- 
mains working only for spatial Wilson loops [16, 17]. This 
parallels percolation properties of monopoles. Moreover, 
it conforms with the observation at high temperatures 
that the spatial Wilson loops keep a string tension in 
contrast to the correlators of Polyakov loops. 

In this paper we will merge the caloron and vortex 
picture focussing on two aspects: (i) to demonstrate how 
the vortex content of iiidiviciual calorons depends on the 
parameters of the caloron solution - in particular the 
holonomy - and (ii) to obtain the vortices in correspond- 
ing caloron ensembles and analyze their percolation prop- 
erties. For simplicity we will restrict ourselves to gauge 
group SU (2) . We will mainly use LCG, which has found a 
correlation of vortices to instantons cores in [18, 19]. We 
recall that LCG has been abandoned for finding vortices 
in SU (2) Monte Carlo configurations because the vortex 
density did not possess a good continuum limit [20] . This 
observation does not invalidate the application of LCG 
to smooth (semiclassical) field configurations. We also 
compare with results obtained by DMCG and IMCG. In 
order to enable the application of these gauge-fixing tech- 
niques we discretize calorons on a lattice, which is known 
to reproduce continuum results very well. 

In the combination of Laplacian Abelian Gauge (LAG) 
and LCG, magnetic monopole worldlines are known to 
reside on the vortex sheets, and we will confirm this for 
the calorons' constituent dyons (see [18] for some first 
findings) . 

In addition, we find another - mainly spatial - part of 
the vortex surfaces. It is strongly related to the twist of 
the dyons within the caloron. We explain this by ana- 
lytic arguments that yield good approximations for the 
locations of these spatial vortices. 

For a single caloron, both parts of the vortex system 
together generate two intersection points needed to con- 
stitute the topological charge in the vortex language.^ 

In caloron ensembles we find the spatial vortex sur- 
faces to percolate only at low temperatures (where the 
holonomy is maximally non-trivial) , while the space-time 
vortex surfaces are rather independent of the phase (i.e. 
the holonomy), both in agreement with physical expec- 



^ Of course, for the caloron as a classical object the topological 
charge density is continuously distributed. 



tations and with observations in caloron gas simulations 
that have evaluated the QQ free energy on one hand and 
the space like Wilson loops on the other [21]. 

From the results for individual calorons it is ck^ar that 
the Polyakov loop, which we treat as an input parameter 
for the caloron solution, is responsible for the percola- 
tion and hence the string tension in the confined phase. 
This lends support for the hypothesis that the holonomy 
is important as the "correct background" for the classi- 
cal objects featuring in a semiclassical understanding of 
finite temperature QCD. 

The paper is organised as follows. In the next two Sec- 
tions II and III we review the properties of calorons and 
vortices, including technicalities of how to discretize the 
former and how to detect the latter. In Section IV we de- 
scribe the vortex content of single calorons. In Section V 
we demonstrate how the vortex content of caloron gases 
changes with the holonomy parameter. We conclude with 
a summary and a brief outlook. Part of our results have 
been published in [22]. 



II. CALORONS 

A. Generalities 

Calorons are instantons, i.e. selfdual"* Yang-Mills fields 
and therefore solutions of the equations of motion, at 
finite temperature. In other words, their base space is 
X where the circle has circumference (3 = l/fc^T 
as usual. 

As it turns out from the explicit solutions [1 3], 
calorons consist of localised lumps of topological charge 
density, which - due to selfduality - are lumps of action 
density, too. For the gauge group SU{N) one can have 
up to N lumps per unit topological charge. When well 
separated, these lumps are static^. Moreover, they pos- 
sess (quantised) magnetic charge equal to their electric 
charge and hence are called dyons. Consequently, the 
moduli of calorons are the spatial locations of the dyons, 
which can take any value, plus phases [10]. 

Another important (superselection) parameter of the 
new solutions by Kraan/van Baal and Lee/Lu [2, 3] is 
the holonomy, the limit of the (untraced) Polyakov loop 
at spatial infinity, 

Voo= lim Vexp(i [ Aodxo] ■ (1) 

\x\^oo \ Jo J 

Due to the magnetic neutrality of the dyons within a 
caloron, this limit is independent of the direction the 



^ The results for antisclfdual calorons with negative topological 

charge are completely analogous. 
^ The gauge field, generically, can and will be time dependent, see 

Sect. II B. 
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limit is taken. (In our convention the gauge fields are 
hermitean, we basically follow the notation of [3] but mul- 
tiply their antihermitean gauge fields by i and reinstate 

/?•) 

In SU{2) we diagonalise Poo, 



Voo = exp (2'!riu!a3) 



(2) 



with (Tj the Pauli matrices. Note that iv = or 1/2 
amount to trivial holonomies "Pqo = ±l25 whereas the 
case oj = 1/4, i.e. trV^o = is referred to as maximal 
nontrivial holonomy. 

The constituent dyons have fractional topological 
charges ("masses") governed by the holonomy, namely 
2u) and 2u) = 1 — 2u!, cf. Fig. 1 upper panel, such that 
- from the point of view of the topological density - the 
constituent dyons are identical in the case of maximal 
nontrivial holonomy u = 1/4. 

To be more concrete, the gauge field of a unit charge 
caloron in the periodic gauge^ is given by 



Vfou^u log 



2770; 



(3) 



/3 



where f] is the 't Hooft tensor (we use the convention 
in [3]) and (f) and % are (a;o-periodic) combinations of 
trigonometric and hyperbolic functions of Xo and x, re- 
spectively, see appendix A and [3]. They are given in 
terms of the distances 

r=\x-zi\, s= \x-Z2\ (4) 

from the following constituent dyon locations 

zt = (0,0,-27rwpV/3), Z2 = (0,0,27rwpV/3) , (5) 

which we have put on the xs-axis with the center of mass 
at the origin (which can always be achieved by space ro- 
tations and translations) and at a distance oi d = irp^ / (i 
to each other. 

In case of large p, the action consists of approximately 
static lumps (of radius P/4ttuj and P/Attlo in spatial di- 
rections) near zi and Z2- In the small p limit the action 
profile approaches a single 4d instanton-like lump at the 
origin. In Ref. [3] one can find more plots of the ac- 
tion density of SU{2) calorons with different sizes and 
holonomies. 

In the far-field limit, away from both dyons the func- 
tion X behaves like 



X 



4d 



|j,g-4wi^r//3g-2-n-ixo _|_ gg-47raJs//3 | 



{r + s + dy 
x[l -|- Q(^Q-^™{'^'^'^^/P'^'^^^/0)y^ ^ 



(6) 




FIG. 1: Action density (top, shown in logarithmic scale and 

cut below e~^'^) and Polyakov loop (bottom) in the (iij.x.j)- 
plane (measured in units of (3) at xq = X2 = for a caloron 
with intermediate holonomy uo — 0.12 and size p — 0.9/3 as 
discretized on a 8 x 48^ x 80 lattice. The dyon locations are 
zi = (0, 0, -0.61) and Z2 = (0, 0, 1.93). 



and hence the off-diagonal part of decays exponen- 
tially, while the Abelian part from 



r + s + d 
r + s — d 



0{e 



-min(47rajr//3,47ra;s / /3) 



(7) 



becomes a dipole field [3] . 

The Polyakov loop in the bulk plays a role similar to 
an exponentiated Higgs field in the gauge group: it is 
+ I2 and — 12 in the vicinity of zi and Z2 ®, respectively, 
cf. Fig. 1 lower panel. The existence of such points is of 
topological origin [24] . Thus the Polyakov loop is a more 
suitable pointer to the constituent dyon locations, which 
agrees with the maxima of topological density for the 
limiting case of well-separated dyons, but is valid even in 
case the two topological lumps merged into one for small 
P- 



^ This gauge is in contrast to the nonperiodic "algebraic gauge" 
where Aq asymptotically vanishes and the holonomy is carried ® On the line connecting the dyons the PolyaJfov loop can actually 
by the transition function. be computed exactly [23]. 



4 



B. The twist 

A less-known feature of the caloron we want to de- 
scribe next is the Taubes twist. It basically means that 
the gauge field^ of one of the dyons is rotated by a time 
dependent gauge transformation (rotated in the direc- 
tion of the holonomy, here the third direction in color 
space) w.r.t. the gauge field of the other dyon when they 
axe combined into a caloron. This is the way the dyons 
generate the unit topological charge [3]. 

The simplest way to reveal the twist is to consider the 
limit of well-separated dyons, i.e. when their distance d 
is much larger than their radius P/Attuj and fS/Airuj. Let 
us consider points near the first dyon, x ~ zi+ 6, where 
the distance S = \S\ is small compared to the separation 
d, but not necessarily compared to the dyon size. Then 
the relevant distances are obviously 

r^6, s= 1(0,0, -d) + (5| = (53 + 0(5 Vd). (8) 

In the appendix we derive the form of the functions 
and X in this limit, 



{x = zi+ S) 



2d 



X{x = zi+S) ~ e 



(5coth(47rcZ)(5//3) - 63 ' 
1 5 



-27rixo / }3 _ 



2rfsinh(47rw(5//3) 



(9) 
(10) 



The large factors of 2d cancel in the expressions log (j) 

and (pd^x relevant for A^, Eq. (3). 

In the vicinity of the other dyon, x = Z2 + 5, with 



s = 5, r = d + 5z + 0{5'^/d) 



(11) 



we get very similar expressions with oj replaced by a; and 
5^ by — (^3, but the time-dependent phase factor is absent, 



(t){x = Z2 + 5) 
X{x = z2 + 5) 



2d 



5coth(47rw5//3) + <53 
1 5 
2rfsinh(47rw(5//5) ' 



(12) 



(13) 



This staticity of course also holds for of this dyon and 
all quantities computed from it. 

Plugging in those functions into the gauge field of 
Eq. (3) one can find that the corresponding gauge field 
components are connected via a PT transformation, and 
the exchange of lo and uj 

(^J*-«^^)(a;o,^2 + 5;w) = -{Al - iAl){-xo,zx-5;i^) 



-4^(a;o,^2 + <?;w) = -Al^{-xo,zi-5;ui) 



(15) 



and a gauge transformation, namely 

A^{xo,Z2 + 6;u)) = -'^Aij,{-xo,zi- S;u!) (16) 
with the time-dependent twist gauge transformation 

(17) 



.Xo 



T{xo) = exp(— TTi— 0-3) . 

This gauge transformation is nonperiodic, T(/3) = — 12 
(but acts in the adjoint representation). 

The Polyakov loop values inside the dyons are obtained 
from x{x = ^1,2 + (5) = 0{S^) and 



(j){x = zi + S) = 
(j){x = Z2 + d) = 



2d 



13/A-Kuj -83 + 0(52) 
2d 

+ 83 + 0((52) 



which results in 
A(i{zi) 

A(i{Z2) 



0-3 







V{z^) = -I2 , 

V{Z2) = +I2 . 



(18) 
(19) 

(20) 
(21) 



Actually, the gauge field around 2*2 is that of a static 
magnetic monopole with the Higgs field $ identified with 
Aq through dimensional reduction. Indeed, it vanishes at 
the core according to (21) and approaches the "vacuum 
expectation value" |$| = 2ttlo/(5 away from the core. Ac- 
cordingly, is identified with DiAq = Fm = Ei, and 
the Bogomolnyi equation with the selfduality equation. 

The gauge field around Zi is that of a twisted 
monopole, i.e. a monopole gauge rotated with T. The 
corresponding Higgs field is obtained from that of a static 
monopole by the same T, transforming in the adjoint rep- 
resentation. Therefore, the Higgs field $ of the twisted 
monopole agrees with the gauge field Aq apart from the 
inhomogeneous term in Eq. (15). $ vanishes at the 
core, too, and approaches the vacuum expectation value 
27rw//3. 

The electric and magnetic charges, as measured in the 
$ direction through the 't Hooft field strength tensor, are 
equal and the same for both dyons. This is consistent 
with the fact that selfdual configurations fulfilling the 
BPS bound must have positive magnetic charge. 

These fields are in some unusual gauge: around the 
dyon cores the Higgs field has the hedgehog form ~ 
{x—Zi^2)a which is called the radial gauge. Far away from 
the dyons the Higgs field i> becomes diagonal up to ex- 
ponentially small corrections. Indeed, if one neglects the 
exponentially small x's of (10) and (13) and replaces the 
hyperbolic cotangent by 1 in the denominator of (9) and 
(12), this would be the so-called unitary gauge with diag- 
onal Higgs field and a Dirac string singularity (along the 
line connecting the dyons). Far away from the caloron's 
dyons, however, the "hedgehog" $ is not "combed" com- 
peletely and there is no need for a singularity^. In other 



The twist can be formulated in a gauge-invariant way by field 

strength correlators between points connected by Schwinger lines ® In contrast, the gauge field A4 written down in Sect. IIA of [9] 
[25]. is diagonal and has a singularity at the xs-axis. 
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words the covering of the color space happens in an ex- 
ponentially small but finite solid angle. 

More precisely, the Higgs field $ approaches 
—2711^0^/(3 and +27ra)(j3//3 away from the static and 
twisting dyon, respectively, for almost all directions. 
These values differ by iray,/ (3, and hence the correspond- 
ing Aqs can be glued together (apart from a gauge singu- 
larity at the origin). Moreover, in Aq the leading far field 
corrections to the asymptotic value, namely monopole 
terms, are of opposite sign w.r.t. the fixed color direction 
(Jz and therefore do not induce a net winding number in 
the asymptotic Polyakov loop. Hence the holonomy is 
independent of the direction. 

We remind the reader that this subsection has been 
dealing with the limit of well-separated dyons, i.e. all 
formulae are correct up to exponential corrections in (3/d 
and algebraic ones in 5/d. 



C. Discretization 

In order perform the necessary gauge transformations 
or diagonalizations of the Laplace operator in numeri- 
cal form we translate the caloron solutions - and later 
caloron gas configurations - into lattice configurations. 
For a space-time grid (with a temporal extent Nq = 8 
and spatial sizes of TV^ = 48, . . . , 80, see specifications 
later) we compute the hnks U^{x) as path-ordered ex- 
ponentials of the gauge field Aij_{x) (for single-caloron 
solutions given by Eqn. (3)). Practically, the integral 

-i J A^{y)dyi,\ (22) 

is decomposed into at least N = 20 subintervals, for 
which the exponential (22) is obtained by exponentiation 
of iAn{y)a/N with A^(y) evaluated in the midpoint of 
the subinterval. These exponential expressions are then 
multiplied in the required order (from x left to x -\- afi 
right). A necessary condition for the validity of this ap- 
proximation is that a/N <^ p with p characterizing the 
caloron size or a typical caloron size in the multicaloron 
configurations. 

Still this might be not sufficient to ensure that the po- 
tential Afj_{y) is reasonably constant within the subinter- 
val of all links and give a converged result. In particular, 
the gauge field (3) is singular at the origin and has big 
gradients near the line connecting the dyons, as visu- 
alised in Fig. 2 of [26]. Hence we dynamically adjust the 
miinber of subintervals for every link, ensuring that 
further increasing N would leave unchanged all entries of 
the resulting link matrix U^{x). 

The lattice field constructed this way is not strictly 
periodic in the three spatial directions, but this is not 
important for the lattices at hand with Ni » A^o- The 
action is already very close to Stt"^, the maximal deviation 
occurs for large calorons {u > 0.9/3) and is about 15 %. 



Lateron, we will make heavy use of the lowest Lapla- 
cian eigenmodes in the LCG. When computing these 
modes in the caloron backgrounds we enforce spatial pe- 
riodicity by hand. In Maximal Center gauges we also 
consider the caloron gauge field as spatially periodic. 



D. Caloron ensembles 

The caloron gas configurations considered later in this 
paper have been created along the lines of Ref. [21]. The 
four-dimensional center of mass locations of the calorons 
are sampled randomly as well as the spatial orientation 
of the "dipole axis" connecting the two dyons and the 
angle of a global C/(l) rotation around the axis in 
color space. The caloron size is sampled from a suitable 
size distribution D(p,T). 

The superposition is performed in the so-called alge- 
braic gauge with the same holonomy parameter oj taken 
for all calorons and anticalorons^. Finally, the additive 
superposition is gauge-rotated into the periodic gauge. 
Then the field Aij,{x) is periodic in Euclidean time and 
possesses the required asymptotic holonomy. We have 
applied cooling to the superpositions in order to ensure 
spatial periodicity of the gauge field. 

In Sect. V we will compare sequences of random 
caloron gas configurations which differ in nothing else 
than the global holonomy parameter to. 



III. CENTER VORTICES 

To detect center vortices, we will mainly use the 
Laplacian Center Gauge (LCG) procedure, which can 
be viewed as a generalization of Direct Maximal Center 
Gauge (DMCG) with the advantage to avoid the Gribov 
problem of the latter [27]. We will compare our results 
to vortices from the maximal center gauges DMCG and 
IMCG in Sect. IV E. 

In LCG one has to compute the two lowest eigenvectors 
of minus the gauge covariant Laplacian operator in the 
adjoint representation^^, 

-A[U^]cPo,i = Ao,i<Ao,i (23) 

+ U^{x - A)'"4-a/i,y - 2^''''<5,,)24) 
a, 6= 1,2,3, 



Superposing (aiiti)calorons with different holonomies would cre- 
ate jumps of Ao in the transition regions. 

We use (j) with a subindex for the eigenmodes of the Laplacian, 
not to be confused with the auxiliary function (j) involved in the 
caloron gauge field, Eq. (3). 



which we do by virtue of the ARPACK package [28]. 

For the vortex detection, the lowest mode (j)o is rotated 
to the third color direction, i.e. diagonalised, 

% = i<^ok3. (25) 

The remaining Abelian freedom of rotations around the 
third axis, V vV with v = cxp^iaa^) is fixed (up to 
center elements) by demanding for <pi a particular form 
with vanishing second component and positive first com- 
ponent, respectively, 

(''^</,i)«=2 = , (''^<^i)«=^ > . (26) 

Defects of this gauge fixing procechire appear when 
(^0 and <j)i are collinear, because then the Abelian free- 
dom parametrised by v remains unfixed. In [27] it was 
shown, that the points x, where 4>o{x) and (f>i{x) arc 
collinear, define the generically two-dimensional vortex 
surface, as the Wilson loops in perpendicular planes take 
center elements. This in chides points x, where (^q van- 
ishes, (j)o{x) = 0, which define monopole worldlines in the 
Laplacian Abelian Gauge (LAG) [29]. 

We detect the center vortices in LCG with the help of 
a topological argument: after having diagonahscd 00 by 
virtue of V, Eq. (25), the question whether cpo and (f>i are 
collinear amounts to ^<j)i being diagonal too, i.e. having 
zero first and second component. We therefore inspect 
each plaquette, take all four corners and consider the 
projections of taken in these points onto the (cti, 0-2)- 
plane, see Fig. 2. 

By assuming continuity-'^ ^ of the field (more pre- 
cisely, its (<Ji, cr2)-projection) between the lattice sites of 
this plaquette, we can easily assign a winding number to 
it. By normalization of the two-dimensional arrows this 
is actually a discretization of a mapping from a circle in 
coordinate space to a circle in color space. In the contin- 
uum this could give rise to any integer winding number, 
while with four discretization points the winding number 
can only take values {—1,0,1}. This winding number 
can easily be computed by adding the angles between 
the two-dimensional vectors on neighbouring sites. 

A nontrivial winding number around the plaquette im- 
plies that the (ci, c72)-components of have a zero 
point inside the plaquette, which in turn means that the 
two eigenvectors are there collinear in color space. In 
this case we can declare the midpoint of that plaquette 
belonging to the vortex surface. The vortex surface is a 
two-dimensional closed surface formed by the plaquettes 
of the dual lattice. The plaquettes of the dual lattice are 
orthogonal to and shifted by a/2 in all directions relative 
to the plaquettes of the original lattice. 

At face value the above procedure is plagued by points 
where the lowest eigenvector is close to the negative 



^ The continuity assumption underlies all attempts to measure 
topological objects on lattices. For semiclassical objects it is 
certainly justified. 
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FIG. 2: The topological argument to detect vortices on a 
given plaquette: The transverse components of the first ox- 
cited mode 4>\ to the direction of the lowest mode 0o (after 
both have been gauge transformed by V) are plotted for the 
four sites of a plaquette. The configuration shown here has 
a nonvanishing winding number, which implies that the two 
eigenvectors are collinear in color space somewhere inside the 
plaquette. 

(T3-direction. Such situations are inevitable when (/)o has 
a hedgehog behavior around one of its zeroes, i.e. for 
monopoles in the LAG. Then the diagonalising gauge 
transformation V changes drastically in space. The cor- 
responding transformed first excited mode may give 
artificial winding numbers and thus unphysical vortices 
if we insist on the continuity assumption in this case. 

Actually, to detect vortices, the lowest eigenvector can 
be fixed to any color direction [27], i.e. to different di- 
rections on different plaquettes. Using this we rotate 
plaquette by plaquette to the direction of the average ^0 
over the four corners of the plaquette. This gauge ro- 
tation is in most cases a small rotation. Afterwards we 
inspect 0i's color components perpendicular to the av- 
erage direction (this can be done by inspecting in 
the ((Ti, (72)-plane after diagonalising the four-site aver- 
aged lowest eigenvector, the resulting gauge transforma- 
tion now changes only mildly throughout the four sites 
of the plaquette) . 

Note that the winding number changes sign under 
<?!)o — > — (^0) but not under — > —<^\ (both changes of 
sign do not change the fact that these fields are eigen- 
modes of the Laplacian). Hence the global signs of (^q^ 
and also the signs of the winding numbers are am- 
biguous. 



IV. VORTICES IN INDIVIDUAL CALORONS 

The following results are obtained for single calorons 
discretized on space-time lattices with Nq = 8 (meaning 
that our lattice spacing is a = /3/8) and Ni = N2 = 
48, N:i = 80 or Ni = N2 = = 64 points. 

For the LCG vortices we have to take an ambiguity into 
account, namely the dependence of the adjoint Laplacian 
spectrum on the lattice discretization, in particular the 
ratio N^/Ni^2- From experience we can summarize that 
the lowest adjoint eigenmode (t)o is rather independent of 
that "aspect ratio" . The first excited mode 01 depends 
on it in the following way, cf. Fig. 3: for large N^/Ni^2 
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FIG. 3: The lowest 30 eigenvalues of the adjoint Laplacian op- 
erator for a caloron with a; = 0.12 and p — 0.7/3 discretized on 
8 X 48^ X 80 (top left) and 8 x 64^ x 64 (bottom left). Two-fold 
degeneracies are plotted as bold lines (and some eigenvalues 
have been slightly shifted to bo distinguishable at this reso- 
lution). For comparison wc plotted in the right panels the 
free spectra on the same lattices marking their degeneracies 
by numbers. The lowest singlets on the rhs. always belong to 
the eigenvalue A = 0. 



the first excited mode is a singlet, whereas for inter- 
mediate and small N^/Ni^2 it is a doublet. 

This ambiguity reflects the fact that we are forcing 

states of a continuous spectrum into a finite volume, 
which ~ like waves in a potential well - are then sen- 



sitive to the periodic boundary conditions"'^^. Localised 
bound states, on the contary, should not depend much 
on the discretization. 

Indeed, the absolute values and degeneracies of the 
eigenvalues can be understood by mimicking the caloron 
with constant links, 

Uo = exp(27ria;c73/iVo) , Ui = l2, (27) 

that reproduce the holonomy (and have zero action). For 
Laplacian modes in the fundamental representation this 
approximation was shown to be useful in [31]. 

In this free-field configuration the eigenmodes are 
waves proportional to H/n 6xp(27rm^a;^/iV^a) with inte- 
ger n^. At nontrivial holonomies and on our lattices with 
Nq <C A^i.2 < jVs one can easily c;onvince oneself, that the 
lowest part of the spectrum is formed by modes in the 
third color direction, ~ 0-3, which do not depend on Xq, 
uq ~ 0. The eigenvalues are then given by trigonomet- 
ric functions of 2TTni/Ni, which for large Ni can be well 
approximated by 

^=^^E(2'^^)' (lowest A). (28) 

In other words, a wave in the ith direction contributes 
nf "quanta" of {2TT/Nif to the eigenvalue. The lowest 

eigenvalue in this approximation is always zero. This fits 
our numerical findings quite well, see Fig. 3. 

In the asymmetric case, = 80, iVi^2 = 48 obvi- 
ously the "cheapest excitation" is a wave along the 0:3- 
axis (connecting the dyons), i.e. n-^ = ±1. This gives 
a doublet, which in the presence of the caloron is split 
into two lines, see Fig. 3 top, the first excited mode is 
thus a singlet (the next modes are those with nontrivial 
ni = ±1 or n2 = ±1 forming an approximate quartet 
and so on). 

In the symmetric case, A^^ = 64, on the other hand, 
excitations along all Xi give equal energy contribution. 
For the excited modes this gives a sextet, which is again 
split by the caloron, see Fig. 3 bottom. It turns out that 
the first excited mode remains two-fold degenerate. The 
eigenmodes are close to combinations of waves with non- 
trivial rii = ±1 and with nontrivial n^. = ±1, reflecting 
the calorons' ax;ial symmetry around the Xs-axis. 

This finally explains the different spectra and different 
shape of the eigenmodes on the different lattices. 



A similar effect lias been observed in Fig. 1 of [30], where the 
adjoint modes in the background of an instanton over the four- 
sphere have been shown to depend on the radius of the sphere. 



8 



A. The lowest eigenvector and the LAG monopoles 

As it turns out, away from the dyons the lowest mode 
(po becomes diagonal^^ and constant, for normahsabil- 
ity reasons it is then approximately (0,0, l/VVol)-^ with 

yoi = N0N1N2N3. 

Near each dyon core we find a zero of the third com- 
ponent of the lowest mode, (f>Q^^, see Fig. 4 top panel. 
Together with the first and second component being very 
small on the whole x^-axis, we expect zeroes in the mod- 
ulus |(/)o| at the constituent dyons, which means that the 
dyons are LAG-monopoles, cf. Fig. 3 in [18] and Fig. 10 
in [31]. 

Such zeroes can be unambiguously detected by a wind- 
ing number on lattice cubes similar to that of Sect. III. 
As a result we find almost static LAG-monopole world- 



lines for large calorons at the locations of their dyons, 
while monopole loops around the caloron center of mass 
are seen for small calorons (with p < 0.5, where the ac- 
tion density is strongly time-dependent as well), see Fig. 
5. Note that these locations are part of the LCG vortex 
surface by definition. Similar monopole worldlines have 
been obtained in the MAG [8, 32]. Adjoint fermionic zero 
modes, on the other hand, detect the constituent dyons 
by maxima [33]. 

The lowest mode also reflects the twist of the caloron: 
the flrst and second component of (f>o near the dyon core 
are either static or rotate once with time Xq evolving 
from to f3. Fig. 6 shows this for the lowest mode 
as well as for the first excited mode. Our results are 
essentially equal to Fig. 9 of [31], just with a resolution 
of A'o = 8 (instead of iVo — 4) more clearly revealing the 
sine- and cosine-like behaviours. 
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FIG. 4: Top: The third component of the lowest mode, 00^^, 
along the a;3-axis (in units of /3 at = /2) for a caloron with 
intermediate holonomy co — 0.1 and size p = 1.0 discretized 
on a 8 X 48^ x 80 lattice. The dyons have xa-locations —0.63/3 
and 2.51/3. Note that for that lattice I/VYoI = 0.00082, a 
value that is indeed taken on by the lowest mode far away 
from the dyons. The other components 0q^^'^ are found to 
be of order 10~* [not shown]. Bottom: the gauge field ^0^^ 
(in units of inverse /3), which is related to the Higgs field $ 
used to explain the behaviour of the lowest mode around the 
dyons (see text). Note that Aq takes the value — 7r//3 near the 
twisting dyon. 



The third direction in color space is distinguished by our (gauge) 
choice of the holonomy, Eqn. (2). 



In order to understand the behaviour found for the 
lowest adjoint mode 4>o, we propose to compare it to the 
Higgs field <I> discussed in Sect. II B. For the static dyon 
one has from time-independence Dq^ = and from the 
equation of motion Di{Di^) = DiFio — 0. Therefore 
<I> of a single static dyon is a zero mode of the adjoint 
Laplacian —A = — D^. For the twisting dyon the same 
equations apply due to the transformation properties of 
<I> (under T) and the latter is again a zero mode of the 
Laplacian. These zero modes approach a constant (the 
vev) asymptotically, so they are normalizable like a plane 
wave. 

Around each dyon core, the lowest adjoint mode 0o be- 
haves similar to $ of that dyon: it vanishes at the dyon 
core, becomes constant and dominated by the third com- 
ponent away from the dyons, it reveals the Taubes twist 
(around the twisting dyon) and is in the same gauge as <&. 
Since the latter are zero modes of — A in the background 
of isolated dyons, a combination of them is a natural 
candidate to be the lowest mode of that (non-negative) 
operator in the caloron background. 

The lowest adjoint mode (j)o for calorons with well- 



. 4 
0.2 



-3 


-2 


-1 


u 


1 


2 


3 


-3 


-2 


-1 


' n ' 



1 


2 


3 



FIG. 5: Zeroes of the lowest adjoint mode, i.e. monopoles 
in Laplacian Abelian Gauge, in the (xo, a;3)-plane (both in 
units of (3, X3 horizontally, at xi = X2 = 0) for calorons of 
holonomy uj — 0.1 and sizes p — 0.5/3 (upper panel, zi = 
(0,0,-0.16), Z2 = (0,0,0.63)) and p = 0.9/3 (lower panel, 
zi = (0,0,-0.51), 22 = (0,0,2.04)). At the origin a closed 
monopole wordline of minimal size occurs, which we ascribe 
to the gauge singularity in the caloron gauge field. 
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FIG. 6: The twist of the caloron gauge field reflected in the 
behaviour of the adjoint Laplacian modes. Shown are the 
individual color components of the lowest mode (top) and 
the first excited mode (bottom) as a function of time xq in 
the vicinity of the twisting dyon at zi. The dashed curve 
depicts the (almost constant) third color component. The 
corresponding plots in the vicinity of the static dyon would 
simply show static lines. 



separated dyons is therefore best described in the fol- 
lowing way, cf. Fig. 4: Around the static dyon at 2*2 
one has (pQ ^ ^ — Aq, where the proportionality con- 
stant of course disappears from the eigenvalue equation 
(23), but is approximately given by the normalization: 
\(j)o\ (0, 0, l/\/Vol)-^. Around the twisting dyon at 
zi, one has to compensate for the inhomogeneous term 
^ = Aq + 0-3 (tt//?) (cf. eqn. (15)). The propor- 
tionahty constant there turns out to be negative, such 
that the lowest mode is able to interpolate between these 
shapes with a rather mild variation throughout the re- 
maining space, see Fig. 4 upper panel. 



B. Dyon charge induced vortex 
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FIG. 7: The dyon charge induced part of the vortex in case 
the first excited mode is a singlet: for a large caloron (top) 
and for a small caloron (bottom), shown schematically in the 
plane xi = X2 = 0. 



plane. In other words, the space-time vortex connects 
the dyons once through the center of mass of the caloron 
and once through the periodic spatial boundary of the 
lattice. 

The magnetic flux (measured through the winding 
number as described in Sect. Ill) at every time slice 
points into the ±a;3-direction. Its sign changes at the 




In the following we present and discuss one part of 
the calorons' vortex that is caused by the magnetic 
charge of constituent dyons. Our findings are summa- 
rized schematically in Figs. 7 and 9. 

The ambiguity of the first excited mode (pi of the ad- 
joint Laplacian influences this part of the vortex most 
such that we have to discuss the singlet and doublet 
cases separately. We find that for the singlet 0i, e.g. 
for N^/Ni^2 — 80/48, the vortex consists of the whole 
(xo, a;3)-plane at a;i = a;2 = only, see Figs. 7 and 8. 
Hence this part of the vortex is space-time like. It in- 
cludes the LAG-monopole worldlines, which are either 
two open (straight) lines or form one closed loop in that 



FIG. 8: The dyon charge induced part of the vortex from 
the singlet first excited mode as measured in a caloron with 
holonomy uj — 0.25 and p = 0.6/3 in a time slice. The outcome 
is identical to the a;3-axis and the same for all time slices. The 
dots denote points on the vortex (lines along xo) where the 
fiux changes, i.e. the LAG-monopoles. 
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FIG. 9: The dyon charge induced part of the vortex from the 
doublet first excited mode for a large caloron (top) and for a 
small (bottom) caloron schematically at X2 ~ 0. 



dyons as indicated by arrows"'^ in Fig. 7. The flux is 
always pointing towards the twisting dyon. 

Independently of the flux one can investigate the 
alignment between the lowest and first excited mode. It 
changes from parallel^^ to antiparallel near the static 
dyon, because the lowest mode 0o vanishes (i.e. the dyon 
is a LAG-monopole) [27]. In addition we find two other 
important facts not mentioned in [27]: the alignment 
does not change at the twisting dyon since both modes 
(/>o and 4>i vanish there and it changes at some other 
locations outside of the calorons' dyons because 0i has 
another zero there [not shown]. 

For the doublet excited mode, i.e. at smaller N-i/Ni,2 = 
64/64, the dyon charge induced vortex is slightly differ- 
ent: again it connects the dyons, but now (for a fixed 
time) via two lines in the "interior" of the caloron, pass- 
ing near the center of mass, see Figs. 9 and 10. These 
lines exist for all times for which the monopole worldline 
exists, that is for all times if the caloron is large and for 
some subinterval of xq if the caloron is small (and the 
monopole worldline is a closed loop existing during the 



We have fixed the ambiguity in the winding number described in 
Sect. Ill by fixing the asymptotic behaviour of the lowest mode. 
In itself, calling 4>0 and (pi parallel is ambiguous as that changes 
when one of these eigenfunctions is multiplied by -1. The tran- 
sition from parallel to antiparallel or vice versa, however, is an 
unambiguous statement. 



subinterval) . 




FIG. 10: The dyon charge induced part of the vortex from 
the doublet first excited state as measured in a caloron with 
holonomy lu — 0.25 and p — 0.6/3 (same as in Fig. 8) at a 
fixed time slice. Like in Fig. 8 the dots denote points on the 
vortex where the flux changes. The a;3-axis has been added 
to guide the eye, it is not part of the vortex surface here. 

These two vortex surfaces spread away from the xa-axis 
which connects the dyons. The axial symmetry around 
this axis is seemingly broken. However, using other linear 
combinations of the doublet in the role of the first excited 
mode (keeping the lowest one) in the procedure of center 
projection, the vortex surface is rotated around the 0:3- 
axis. The situation is very similar to the "breaking" of 
spherical symmetry in the hydrogen atom by choosing a 
state of particular quantum number m out of a multiplet 
with fixed angular momentum /. 

The magnetic fiux fiips at the dyons, just like in the 
case with singlet 

Notice that these vortices are predominantly space- 
time like, but have parts that are purely spatial, in par- 
ticular for small calorons, namely at minimal and maxi- 
mal Xq of the dyon charge induced vortex surface (and at 
other locations in addition, when the smooth continuum 
surface is approximated by plaquettes). 

C. Twist-induced vortex 

In this section we will discuss the second part of the 
LCG vortex surfaces we found for individual calorons. 
We start again by discussing the singlet case. The twist- 
induced vortex in the singlet case appears at a fixed time 
slice and hence is a purely spatial vortex. In contrast to 
the space-time part, this vortex surface does not contain 
the monopole/dyon worldlines. Hence it is not obvious 
that this part of the vortex structure is caused by them. 
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FIG. 11: Twist-induced part of the vortex ("bubble") from singlet first excited modes for calorons of size p = 0.6/9 and 
holonomies from left to right: uo — 0.1,0.12,0.16 (upper row) ui = 0.2,0.25,0.3 (middle row) and cj — 0.34 (lower row, left 
panel). The plot in the lower right panel summarises the results for uj — 0.1, 0.12, 0.16, 0.2, 0.25, 0.3, 0.34, at xi — 0, i.e. the 
bubbles are cut to circles. The plane near the boundary in the ui = 0.25 picture is an artifact caused by periodic boundary 
conditions. 



The properties of this spatial vortex depend strongly 
on the holonomy, which will be very important for the 
percolation of vortices in caloron ensembles in Sect. V. 

In short, our finding is that the twist-induced part of 
the vortex is a closed surface around the twisting dyon as 
long as the holonomy parameter is w < 1/4, and becomes 
a closed surface around the static dyon for uj > 1/4, we 
will refer to these surfaces as "bubbles" . For maximal 
nontrivial holonomy uj — 1/4: the vortex is the = 
plane, i.e. the midplane perpendicular to the axis con- 
necting the dyons, we will refer to it as "degenerate bub- 



ble" . 

The bubble depends on the holonomy lo as shown in 
Fig. 11. For two complementary holonomies uj = luq and 
u! = ^ — ujQ the bubbles are of same shape just reflected 
at the origin, thus one of them encloses the static dyon 
and another encloses the twisting dyon. This is to be 
expected from the symmetry of the underlying calorons. 
In the limit of w ^ 1/4 the bubbles grow to become a 
flat plane which enables to turn over to the other dyon. 

In our UJ = 0.25 data we find another piece of the vortex 
near the boundary of the lattice, see Fig. 11. It is an 
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FIG. 12: Spatial part of the vortex ("bubble") for calorons of fixed intermediate holonomy = 0.12 and sizes from left to right: 
p = 0.6/3, 0.7/3, 0.9/3. The panel on the very right shows a summary of the bubbles for p = 0.6/3, 0.7/3, 0.8/3, 0.9/3 at a;i = 0. 
That the bubble for p = 0.9/3 is much bigger than that for p — 0.7/3, 0.8/3 is probably a finite volume effect. For small sizes 
p (and also in the limiting cases of holonomy lj close to the trivial values and 1/2) we have met difficulties in resolving the 
corresponding small bubbles in the lattice dicretization. 



artefact of the finite periodic volume. Likewise, very large 
bubbles in our results have deformations since they come 
close to the boundary of the lattice. The intermediate 
bubbles shown in these figures are generally free from 
discretization artefacts and can easily be extrapolated 
(at least qualitatively) to these limits. 

The size of the bubble also depends on the size param- 
eter p of the caloron, i.e. the distance between the dyons, 
as shown in Fig. 12. 

The time-coordinate of LCG bubbles in large calorons 
is always consistent with xq — 0.5/3. For small calorons, 
on the other hand, xq = is the exclusive time slice: the 
action density peaks there and the LAG monopoles are 
circling around it (cf. Fig. 5). However, the bubbles of 
small calorons are too small to be detected. 

In the case of the first excited mode being a doublet, 
similar bubbles have been found. They also enclose one 
of the dyons and degenerate to the midplane for to — 
1/4. Their sizes, however, may be different and they 
are distributed over several time slices. Considering the 
collection of all time slices, these fragments add up to full 
bubbles. 



1. Analytic considerations 

In the following we present two analytic arguments - 
relying on the twist - that support the existence of the 
bubbles (playing the role of spatial vortices) and help to 
estimate their sizes. 

The first one is specific for vortices in LCG. As we have 
demonstrated in Sect. IV, the lowest mode (t)Q twists near 
the twisting dyon and is static near the static dyon. The 
same holds for the first excited mode (j>i, see Fig. 6. 

Then a topological argument shows that they have to 
be (anti)parallel somewhere inbetween, cf. Fig. 13. As 
(j)o, 4>i and the diagonalising gauge transformation V 



are static around the static dyon, so is (pi and its pro- 
jection along the third direction (see the right part of 
Fig. 13). We assume that this projection is nonzero, oth- 
erwise the two states are obviously (anti)parallel and the 
point would belong to the vortex already^®. 

In the twisting region called S, the two lowest modes 




FIG. 13: Behaviour of the nondiagonal elements of V tran- 
formed first excited mode ^(f>i in the twisting region (left) and 
in the static region (right) with time xo evolving upwards. 
The lattice sites inbetween are indicated only at xq = 0. On 
the entire discretised rectangle the field has winding number 1, 
meaning it contains the twist-induced vortex. More precisely, 
it is the "plaquette" marked with filled circles that contains 
the winding (in analogy to Fig. 2) and thus the vortex (in all 
other plaquettes the field performs a partial winding but then 
winds back). 



in particular to the space-time part since then the two modes are 
(anti)paraUel for all xq 
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behave like (suppressing arguments x) 



static region. 



(poixa) = T{xo)(I)q{xo = 0)T^(a;o) 
(piixa) = T(a;o)(/)i(xo = 0)T^(a;o) 



(29) 



with the twisting transformation/rotation from Eq. (17). 
The time dependence of the diagonahsing V can be de- 
duced easily^'', 



such that 



F(a;o) = T(xo)nO)T^(a;o) 



^0i(xo)=T(xo)%(O)Tt(a;o) 



(30) 



(31) 



Again we assume that the two modes are not 
(anti) parallel at xq = 0. Then ^(/)i(0) has a nonvanishing 
component perpendicular to CT3. According to Eq. (31) 
this component then rotates in time Xq around the third 
direction (left part of Fig. 13). This immediately implies 
that there is a space-time "plaquette" (in the sense of 
Fig. 2, marked in Fig. 13 with filled circles) that contains 
a point where the two modes are coUinear. Notice the 
similarity of Figs. 13 and 2. 

This argument applies to all pairs of points with one 
point in the twisting region S and one point in the static 
region (its complement) S: on any line connecting the 
two there exists a point which belongs to the vortex. This 
results in a closed surface at the boundary between S and 
S (see below) . The time-coordinate of this surface is not 
determined by these considerations. 

Our argument can be extended to vortices beyond 
LCG. For that aim we mimic the caloron gauge field 
by = 0, Ai = in the static region S and Aq = 
— cr3(7r//3), = in the twisting region S (cf. Eq. (15)) 
[22] . In this simplified gauge field vortices can be located 
directly by the definition that — 1 Wilson loops are linked 
with them. Obviously rectangular Wilson loops connect- 
ing (0,fi), (/3,fi), {/3,X2), (0,^2) and (0,fi) are_-l2 
if and only if belongs to S and X2 belongs to S (or 
vice versa). This again predicts spatial vortices at the 
boundary between the twisting and the static region. 

Actually, this argument is exact if one chooses for the 
points xi_2 the dyon locations 21^2: the path ordered ex- 
ponentials at fixed xi^2 are the Polyakov loops =pl2 and 
the remaining spatial parts are inverse to each other be- 
cause of periodicity and cancel. Hence there should al- 
ways be a spatial vortex between the two dyons. 

Thus the twist in the gauge field of the caloron itself 
gives rise to a spatial vortex. This vortex extends in the 
two spatial directions perpendicular to lines connecting 
S and S, just like a bubble. 

Note that the two arguments above do not work 
purely within the twisting region or purely within the 



^"^ The first factor is necessary, otherwise V is singular around the 
north pole and nonperiodic. 



It remains to be specified where the boundary between 
the twisting region S and its complement S is. To that 
end one should consider the competing terms - twisting 
vs. static - in the relevant function x, see Eq. (A. 3). Ac- 
tually its derivatives enter the off-diagonal gauge fields, 
see Eq. (3). In the periodic gauge we have used so far, 
there is an additional term proportional to x itself. To 
decide whether the static or the twisting part dominates 
(at a given point) it is better to go over to the algebraic 
gauge, where this term is absent and where x must be re- 
placed by X = exp{4,Triujxo/ l3)x [3]. The two competing 
terms become 



„4^^u^x„/0 Smh{4TTUjr/(3) _ 

tpr 

-47r»i;xo/;9 sinh(47UJs//3) _ 



fst 



ft\ 



(32) 
(33) 



with -0 given in Eqn. (A.l). Note that the time de- 
pendence of these functions still differs by a factor 
exp(27rza:o//3). 

We finally define the twisting region S as where the 
gradient of /twist dominates 



/J. ^ static I 



(34) 








FIG. 14; The bubbles measured for calorons with holonomy 
u) = 0.12, size p — 0.6/3 (top) and p = 0.9/3 (bottom) respec- 
tively as a function of X2 (vertically) and X3 (horizontally) at 
xi = compared to the boundary of the twisting region S, 
the smooth curve computed from the equality in Eq. (34). 
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and its boundary where the equahty holds. 

In two particular cases this can be determined ana- 
lytically. For the case of maximally nontrivial holonomy 
uj = uj = 1/4, the two functions /static. twist only differ 
by the arguments r vs. s. Then the boundary of S is 
obviously r = s, which gives the midplane between the 
dyons. This indeed amounts to our numerical finding, 
the degenerate bubble for ui = 1/4. 

In the large caloron limit and if we further as- 
sume the solutions of the equality in Eq. (34) to obey 
ujr/(3,ujs//3 > 1, it is enough to compare in /static,twist 
the exponentially large terms in sinh and -0. This yields 
for the boundary of S the equation ur = lus, which can 
be worked out to give 

xj + xl + ixa-ndf = (nd)^ , n = ^^ (35) 

CO — CO 

Thus, the boundary of the twisting region S is a sphere 
with midpoint (0, 0, ild) and radius \il\d. This sphere al- 
ways touches the origin, is centered at negative and pos- 
itive X3 for w < 1/4 and > 1/4, respectively, and again 
degenerates to the midplane of the dyons for maximally 
nontrivial holonomy lo ~ Cu — 1/4. 

In Fig. 14 wc compare the boundary of S obtained 
from the equality in Eq. (34) to the numerically obtained 
bubbles in LCG for two diS^erent values of the caloron 
parameter p. The graphs agree qualitatively. 

One could also think of characterising the locations x 
of the twist-induced vortex by a fixed value of the traced 
Polyakov loop, say trV{x) = 0. This also encloses one 
of the dyons and becomes the midplane for w = 1 /4. In 
the large separation limit, however, this surface is that 
of a single dyon of fixed size set by /3 and lu (just like the 
topological density) . It does not grow with the separation 
d, which however seems to be the case for the measured 
vortices as well as for the boundary of S using the far field 
limit, Eq. (35). Hence the local Polyakov loop seems not 
a perfect pointer to the spatial vortex. 



D. Intersection and topological chsirge 

To a good approximation the dyon charge induced vor- 
tex extends in space and time connecting the dyons twice, 
whereas the twist induced vortex is purely spatial around 
one of the dyons. This results in two intersection points 
generating topological charge as we will describe now. 

The notion of topological charge also exists for (sin- 
gular) vortex sheets. In order to illustrate that let us 
choose a local coordinate system and denote the two di- 
rections perpendicular to the vortex sheet, in which a 
Wilson loop is —1, by fi and u. The Wilson loop can 
be generated by a circular Abelian gauge field decaying 
with the inverse distance, which generates a gauge field 
F^i, (the magnetic field, say B3 oc F12 for a static vortex 
in the a;3-direction, is tangential to the vortex, respec- 
tively). The corresponding fiux is via an Abelian Stokes' 



Theorem connected to the Wilson loop and is nothing but 
the winding number used in LCG to detect the vortex. 

In order to generate topological charge proportional 
to e^"'"' F^^Fp^^, the vortex thus needs to "extend in all 
directions". This is made more precise by the geomet- 
ric objects called writhe and self-intersection. The re- 
lation to the topological charge including example con- 
figurations has been worked out for vortices consisting 
of hypercubes in [34] and for smooth vortices in [35]. 
The result is that a (self) intersection point - where two 
branches of the vortex meet such that the combined tan- 
gential space is four-dimensional - contributes ±1/2 to 
the topological charge. The contribution of the writhe 
is related to gradients of the vortex' tangential and nor- 
mal space w.r.t. the two coordinates parametrising the 
vortex. Two trivial examples are important for vortices 
in a caloron: a two-dimensional plane as well as a two- 
dimensional sphere embedded in four-dimensional space 
have no writhe. Since the two parts of our vortex are of 
these topologies, we immediately conclude that the topo- 
logical charge of vortices in calorons comes exclusively 
from intersection points. 

We first discuss the position of the intersection points 
in the singlet case, cf. Fig. 15 top panel. The twist in- 
duced bubble occurs at a fixed time slice and so does any 
intersection point. The dyon induced vortex consists of 
two static straight lines from one of the dyon to another 
and therefore intersects the bubble twice on the .I'u-axis. 

There are two exceptions to this fact: for small 
calorons (small p) the dyon induced vortex exists in some 
time slices only and the number of intersection points 
depends on whether the time-coordinate of the bubble is 




FIG. 15: The intersection of the spatial bubble (at fixed xa = 
13/2) for < 1/4 witli the space-time part of the vortex from 
the singlet (top) and the doublet (bottom) excited mode (in 
the doublet case the bubble is distributed over several time 
slices). 
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within that timc-intcrval (but the bubble is usually too 
small to detect when p is small). 

The case of maximal nontrivial holonomy is particular 
because for the corresponding degenerate bubble there is 
only one intersection point at the center of mass of the 
caloron, the other one moved to 0:3 — > ±00 as a; — > l/4±0 
in infinite volume. 

Concerning the sign of the contributions, it is essential 
that the relative sign of the vortex flux is determined: 
the magnetic flux of the dyon induced vortex flips at 
the dyons and hence is of opposite sign at the intersec- 
tion points on the bubble. One can depict the flux on 
the bubble by an electric field normal to the bubble (i.e. 
hedgehog-like). It follows that in LCG the contributions 
of the intersection points to the topological charge of the 
vortex are both +1/2. 

The vortex in the caloron is thus an example for a 
general statement, that a non-orientable vortex surface 
is needed for a nonvanishing total topological charge. In 
our case the two branches of the dyon charge induced 
vortex have been glued together at the dyons in a non- 
orientable way: the magnetic fluxes start or end at the 
dyons as LAG-monopoles (this construction is impossible 
for the bubble as the dyons are not located on them). 
Thus, vortices without monopoles on them possess trivial 
total topological charge. 

In the doublet case with its fragmented bubbles there 
axe still two intersection points (cf. Fig. 15 bottom panel) 
which again contribute topological charges of +1/2 each. 

To sum up this section we have demonstrated that the 
vortex has unit topological charge like the caloron back- 
ground gauge field. This result is not completely triv- 
ial as there is to our knowledge no general proof that 
the topological charge from the gauge field persists for 
its vortex "skeleton" after center projection (P- vortices) . 
Moreover, the topological density of the caloron is not 
maximal at the two points where the topological density 
of the vortex is concentrated and the total topological 
charge of the caloron is split into fractions of 2uj and 2uj 
whereas that of the vortex always comes in equal frac- 
tions 1/2 from two intersection points, close to the static 
dyon if a; < 1/4 and close to the twisting dyon if w > 1 /4. 

E. Results from Maximal Center Gauges 

We have performed complementary studies of vortices 
both in the Direct and in the Indirect Maximal Center 
Gauges (DMCG [36] and IMCG [37], respectively). The 
DMCG in SU{2) is defined by the maximization of the 
functional 

Ft,mcg[U] = (tr <'U^{x)f , (36) 

with respect to gauge transformations g{x) e SU{2). 
U^{x) are the lattice links and ^U^{x) = g^{x)U^{x)g{x+ 

fj) the gauge transformed ones. Maximization of (36) 
minimizes the distance to center elements and fixes the 



gauge up to a Z{2) gauge transformations. The corre- 
sponding, projected Z{2) links are defined as 

Z^(a;) = sign(tr9[/^(x)) . (37) 

The Gribov copy problem is known to spoil gauges with 
maximizations such as DMCG. In practice we also ap- 
plied random SU{2) gauge transformations before max- 
imizing Fdmcg and selected the configuration with the 
largest value of that functional. 

The IMCG goes an indirect way. At first, one fixes the 
maximal Abelian gauge (MAG [38]) by maximizing the 
functional 

Fmag[U] = ^tr {'^U^{x)a3{sU^{x))^a3) , (38) 

with respect to gauge transformations g € SU{2). The 
MAG minimizes the off-diagonal elements of the links 
and fixes the gauge up to U{1). Therefore, the following 
projection to a C/(l) gauge field through the phase of the 
diagonal elements of the links, 9^{x) = arg ((^i7^(a:;))^-^), 
is not unique. Exploiting the remaining U{1) gauge 
freedom, which amounts to a shift 9fj,{x) — > "On{x) = 
—a{x) + 6n{x) + a{x + fi), one maximize the IMCG func- 
tional 

FimcgIU] = ^ {cos{"e^{x))f , (39) 

IJ.,X 

that serves the same purpose as Fdmcg in (36). Finally, 
the projected Z{2) gauge links are defined as 

Z^(x) = sign(cos("^^(a;))) . (40) 

Finally, th(^ Z{2) links are used to form Z{2) plaque- 
ttes. All dual plaquettes of the negative Z{2) plaquettes 
form closed two dimensional surfaces - the vortex 
surfaces. 

In the background of caloroiis wc tried to confirm both 
dyon charge induced and twist-induced vortices seen in 
LCG. In DMCG, the dyon charge induced vortices are 
observed and the twist induced part splits into several 
parts in adjacent time slices. Choosing the best among 
random gauge copies, the dyon charge induced part dis- 
appears and the twist-induced vortex bubble occurs at 
fixed time slice. In both cases, the bubble is much smaller 
than that found in LCG. 

In IMCG, the situation for dyon charge induced vor- 
tices is rather stable: we find always a space-time vor- 
tex connecting the dyons or better to say the Abelian 
monopoles representing the dyons in MAG [8, 32]. Sim- 
ilar to the LCG doublet case two vortex lines pass near 
the center of mass of the caloron. This structure propa- 
gates either statically or nonstatically in time, depending 
on the distance between dyons in the caloron. 

Straight lines through the center of mass and through 
the outer space, as found in the LCG singlet case, were 
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never observed. One can convince oneself, that it is ac- 
tually impossible to get such vortex structures from Z(2) 
link configurations. 

The situation with twist-induced vortices is unsta- 
ble, as a rule they do not appear in IMCG. The rea- 
son for this could be partially understood in IMCG 
considerations. Let us consider two gauge equivalent 
Abelian configurations that generate local Polyakov loops 
^tr'P{x) = cos(a(.T)) and how their temporal links con- 
tribute to the Abelian gauge functional i^iMCG given in 
Eq. (39). In the quasi-temporal gauge when all sub- 
sequent temporal links are the same, they give a con- 
tribution equal to iVo cos(a(x)/iVo)^ to the correspond- 
ing part of the functional. When, on the other hand, 
all but one temporal links are trivial, the contribution 
is equal to A^o — 1 + cos(a(x))^. For cos{a{x)) ^ —1 
and for sufficiently large A'o we have Aq cos(o(.x)/A"o)^ ^ 
A^o - a{xY/No > No - sin(o(f))^. This means that af- 
ter maximizing the functional (39) and projecting onto 
Z{2), we get all temporal links trivial in all points x 
where the Polyakov loop is not equal to —1. So, the 
twist-induced vortex shrinks to one point where the (un- 
traced) Polyakov loop is equal to — 12. 

Maximization of the functional (39) is equivalent to 
the minimization of the functional 

i^[C/] = ^(sin(«^4a;)))^ (41) 
If we would replace it by the functional 

the situation with the Z{2) projected Polyakov loop 

would change because A^o i/sin(a(.x)/A^o)^ > \/sin(a(x))2 
on points where cos(a(.T)) < and now we would have 
— 12 temporal links in some time slice and trivial tem- 
poral links in all other time slices in the spatial region 
where the Polyakov loop is negative as well as trivial 
temporal Z{2) links in all time slices in the region where 
the Polyakov loop is positive. Numerical studies support 
the appearance of twist-induced vortex on the boundary 
where initial Polyakov loop changes the sign from nega- 
tive to positive. 

One may conclude from this section that the Gri- 
bov copy problems of DMCG and IMCG persist for the 
smooth caloron backgrounds. The basic features of the 
vortices can be reproduced, but for clarity we stick to the 
vortices obtained in the Laplacian Center gauge. 

V. VORTICES IN CALORON ENSEMBLES 

In this section we present the vortex content of ensem- 
bles of calorons. The generation of the latter has been 
described in Sect. II D. We superposed 6 calorons and 6 
anticalorons with an average size of p = 0.6/3 on a 8 x 64^ 
lattice. 



The most important feature of these ensembles is their 
holonomy Vao = exp (27ria;(73). Under the conjecture 
mentioned in the introduction we will use ^tvVoo = 
cos(27ra;) as equivalent to the order parameter {^trV). 
In particular, caloron ensembles with maximally nontriv- 
ial holonomy u = 1/4 mimic the confined phase with 
(itrP) = 0. 

For each of the holonomy parameters lo = 
{0.0625, 0.0125, 0.01875, 0.25} we considered one 
caloron ensemble with otherwise equal parameters. 
Again we computed the lowest adjoint modes in these 
backgrounds and used the routines based on winding 
numbers to detect the LCG vortex content. 

In Figs. 16 and 17 we show the space-time part respec- 
tively the purely spatial part of the corresponding vor- 
tices as a function of the holonomy oj. One can clearly see 
that with holonomy approaching the confining value 1/4 
the vortices grow in size, especially the spatial vortices 
start to percolate, which will be quantified below. 

Fig. 18 shows another view on this property. In these 
plots we have fixed one of the spatial coordinates to 
a particular value, such that vortices become line-like 
or remain surfaces (and may appear to be non-closed, 
when they actually close through other slices than the 
fixed one). These plots should be compared to Fig. 7 
of Ref. [17], which however does not show a particular 
vortex configuration, but the authors' interpretation of 
measurements (in addition, the authors of [17] seem to 
have overlooked that vortices cut at fixed spatial coordi- 
nate still have surface-like parts, i.e. dual plaquettes). 

We find that vortices in the deconfined phase tend 
to align in the time-like direction, while in the confined 
phase vortices percolate in the spatial directions. We re- 
mind the reader that we distinguish the different phases 
by the values of the holonomy, uj ~ 0, w ~ 1/2 vs. 
OJ ^ 1/4, and not by different temperatures, which would 
lead to different caloron density and size distribution. 

We start our interpretation of these results by the fact 
that the caloron background is dilute in the sense that the 
topological density is well approximated by the sum over 
the constituent dyons of individual calorons (of course, 
the long-range A"^^ components still "interact" with the 
short-range A'^^^''^ components inside other dyon cores). 
Therefore it is permissible and helpful to interpret the 
vortices in caloron ensembles as approximate recombi- 
nation of vortices from individual calorons presented in 
the previous sections. 

Indeed, the space-time vortices resemble the dyon- 
induced vortices which are mostly space-time like and 
therefore line-like at fixed xq- The spatial vortices, on 
the other hand, resemble the twist-induced bubbles in 
individual calorons. 

Following the recombination interpretation, the bub- 
bles should become larger and larger when the holon- 
omy approaches the confining holonomy to = 1/4, where 
they degenerate to flat planes. This is indeed the case in 
caloron ensembles: towards cij = 1/4 the individual vor- 
tices merge to from one big vortex, see also the schematic 
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FIG. 16: The space-time part of vortices in caloron ensembles in a fixed time slice. Only the holonomy varies from left to right; 
uj — {0.0625, 0.0125, 0.01875, 0.25} (from deconfined phase to confined phase). The upper row shows the entire vortex content 
in each caloron ensemble, the lower row shows the corresponding biggest vortex cluster. 




FIG. 17: The spatial part of vortices in the caloron ensembles of Fig. 16 (with the same values of the holonomy uj) summed 
over all time slices. Again the upper row shows the entire vortex content and the lower row the corresponding biggest vortex 
cluster. 



plot Fig. 7 in [22]. In other words, the maximal nontrivial 
holonomy has the effect of forcing the spatiall vortices to 
percolate. Consequently the vortices yield exponentially 
decaying Polyakov loop correlators, the equivalent of the 
Wilson loop area law at finite temperature (see below) . 

Note also that there is no similar scenario for the 



space-time vortices. The dyon-induced vortices in each 
caloron are either always as large as the lattice (in the 
singlet case) or are always confined to the interior of the 
caloron (in the doublet case). This is consistent with 
the physical picture, that spatial Wilson loops do not 
change much across the phase transition. 
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holonomy parameter 0.0625 0.125 0.1875 0.25 



space-time extension 
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56 


56 56 


spatial extension 
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56 56 
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FIG. 18: The vortex content of a caloron ensemble in the 
deconfined phase (upper panel, mimicked by holonomy u: = 
0.0625 close to trivial) and in the confined phase (lower panel, 
maximal nontrivial holonomy u = 0.25) in a lattice slice at 
fixed spatial coordinate. The short direction is xq while the 
other directions are the remaining spatial ones, all given in 
units of (5. 



In order to quantify the percolation we measured two 
observables. The first one concerns the spatial and space- 
time extensions of the largest vortices. We have included 
plots of them in the second row of Figs. 16 and 17, re- 
spectively. 

For the spatial extension we superpose the purely spa- 
tial vortex plaquettes of all time slices in one 3d lattice, 
whereas for the space-time extension we remove all purely 
spatial vortex plaquettes. Then we pick the largest con- 
nected cluster in the remaining vortex structure. 

Table I shows the extension of largest spatial and 
space-time vortex clusters for different holonomy pa- 
rameters of otherwise identical caloron ensembles. The 
spatial-spatial vortex cluster extension changes drasti- 



TABLE I: Extensions of the largest vortex cluster (see text) in 
the caloron ensembles of Figs. 16 and 17. Note that the largest 
extension on a 8 • 64^ lattice is \/(8/2)2 + 3- (32/2)2 ^ 55 5 



cally with the holonomy parameter a; whereas the space- 
time one almost keeps to percolate. 

The second row in Table I is related to confinement 
generated by vortices. If center vortices penetrate a Wil- 
son loop with extensions T and L randomly, the proba- 
bility to find n vortices penetrating the area A = TL is 
given by the Poisson distribution 

Pr{n;T,L)=^^^e-P^ (43) 
n! 

where p is the density of (spatial) vortices. The index 
r characterizes the perfect randomness of this distribu- 
tion, in order to distinguish it from an arbitrary empirical 
distribution. The average Wilson loop in such a center 
configuration is given by an alternating sum 

(W^(r,L))=^(-l)"P(n;r,i) (44) 

n 

Obviously, for the Poisson distribution Pr one obtains an 
area law, log{W) oc TL, with a string tension a = 2p. 

We explore space-time Wilson loops {W{L,P)), which 
amount to Polyakov loop correlators at distance L and 
probe confinement, as well as purely spatial Wislon loops 
as a function of their area, {W{L, L')) = {W{A = LL')). 
As Fig. 19 clearly shows, the values of log(M^(L, j3)) show 
a confining linear behaviour (like from random vortices) 
reaching larger and larger distances L when the holonomy 
approaches the confinement value oj — 1/A. At the same 
time the corresponding string tensions also grow by an 
order of magnitude. 

One can also explain the deviation from the confining 
behaviour in holonomies far from 1/4 . The probability 
P(2; L, (3) for two vortices penetrating the Wilson loop is 
found much larger than in the case of Poisson distribution 
[not shown]. This comes from small bubbles, which very 
likely penetrate a given rectangular twice. 

The behaviour of log(W^(A)) on the other hand changes 
only slightly for different holonomies, see Fig. 20. The 
corresponding slopes ("string tensions") vary by a factor 
of approximately 2. For holonomy uj — 0.25 we find an 
exponential decay stronger than proportional to the area. 
A quantitative analysis of this effect needs to include suit- 
able caloron densities and size distributions around the 
critical temperature. 
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FIG. 19; The observable — log{W{L, (5)) as a function of L in units of lattice spacing a = /3/8 from vortices in caloron ensembles 
with holonomies as in Figs. 16 and 17 (from left to right: uj — {0.0625, 0.0125, 0.01875, 0.25}). Note the very different scales. 
Percolating vortices, i.e. random penetrations, would give a linear behavior. The line shown here is a linear extrapolation of 
the first data point. The caloron vortices follow this line for higher and higher distances L when the holonomy approaches the 
confinement value 1/4, thus inducing confinement in the Polyakov loop correlator (see text). 




FIG. 20: The observable ~log{W{A)) as a function of the 
area A in lattice units — I3'^/6A from vortices in caloron 
ensembles with holonomy parameters uj — 0.0625 (left) and 
UJ = 0.25 (right). Percolating vortices, i.e. random penetra- 
tions, would give a linear behavior. The line shown here is a 
linear extrapolation of the first data point. 



VI. SUMMARY AND OUTLOOK 

In this paper we have extracted the vortex content 
of SU(2) calorons and ensembles made of them, mainly 
with the help of the Laplacian Center Gauge, and 
studied the properties of the emerging vortices. Our 
main results are 

(1) The constituent dyons of calorons induce zeroes 
of the lowest adjoint mode and therefore appear as 
monopoles in the Laplacian Abelian Gauge. The cor- 
responding worldlines are either two static lines or one 
closed loop for large and small calorons, respectively. 

(2) One part of the caloron's vortex surface contains 
the dyon/monopole worldlines. The vortex changes 
its flux there (hence the surface should be viewed as 
non-orientable) . These are general properties of LCG 
vortices. The specific shapes of these dyon-induced 
vortices depend on the caloron size as well as on the 
lattice extensions. These vortices are predominantly 
space-time like. 

(3) Another part of the vortex surface consists of a 
"bubble" around one of the dyons, depending on the 
holonomy. The bubble degenerates into the midplane of 
the dyon "molecule" in the case of maximal nontrivial 
holonomy uj = 1/4. This part is predominantly spatial. 
We have argued that it is induced by the relative twist 



between different dyons in the caloron. 

(4) Both parts of the vortex together reproduce the 
unit topological charge of the caloron by 2 intersection 
points with contributions 1/2. 

(5) In dilute caloron ensembles ~ that differ only in 
the holonomy mimicking confined and deconfined phase 
- the vortices can be described to a good approximation 
by recombination of vortices from individual calorons. 
The spatial vortices in ensembles with deconfining 
holonomies a; ~ form small bubbles. With the 
holonomy approaching the confinement value lu — 1/4, 
the spatial bubbles grow and merge with each other, i.e. 
they percolate in spatial directions. We have quantified 
this by the extension of the largest cluster on one hand 
and by the quark-antiquark potential revealed by the 
Polyakov loop correlator on the other. 

In particular the last finding is in agreement with the 
(de)confinement mechanism based on the percolation of 
center vortices. 

We postpone a more detailed analysis of this mech- 
anism as well as working out the picture for the more 
physical case of gauge group SU{3) to a later paper. 
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Appendix 
A. Caloron gauge fields 

Here we give the functions necessary for the gauge 
fields of calorons [3] including (3 and their form in the 
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limits described in Sect. II. The first set of two auxiliary dimensionless functions are 



J.2 _|_ g2 _j_ ^2 

^ = — cos(27r.xo//3) + cosh (47ra;r//3) cosh (47ra;s//3) H sinh(47ra)r//?) sinh(47rws//3) 

/sinh(47rwr//3) , , , ^, sinh(47ra;s//3) \ ,. 
+d ( ^— — ^ cosh(47ra;s//3) + cosh(47ra;r//3) — ^ j (A.l) 

7.2 _|„ g2 _ ^2 

^ = - cos(27rxo//3) + cosh(47ra;r//3) cosh(47ra;s//3) H sinh(47rwr//3) sinh(47rws//3) . (A. 2) 

We remind the reader that r = \x—zi\ and s = \x—Z2\ are the distances to the dyon locations and d = — Z2I — t^P^ / P 
is the distance between the dyon locations, the "size of the caloron" . The next set of two auxiliary functions entering 
Eqn. (3) are 

V' ~, _ 1 ^ /^sinh(47ra;r//3) , ^_2^i^„/^sinh(47ra;s//3) 



For the twist we have analyzed the limit of large size 
d » /? in Sect. II. For points x = zi + S near the location 
of the first dyon, r = \S\ is small and s = d — S3 + 0{S'^ /d) 
is large. Hence the argument Attojs/P is much larger than 
1 (unless trivial holonomy u) = 0) and the hyperbolic 
functions can be replaced by exponential functions with 
exponentially small corrections. On the other hand, no 
manipulations are made in all functions with argument 
4:7ru)r/(3, such that we get the exact expressions in terms 
of the distance r, 



i> = \^e'^-^'^'l<^ ('cosh(47ra)r//3) - - sinh(47ra;r//?) 



d 



(A.4) 



The exponentially large prefactors cancel in the functions 

and X- 



(j){x = Zi + S) 
Xix = zi + 5) 



2d 

5coth(47ra)^//3) - ^3 

_ .-27ria;o//3i 1_ 



2dsinh(47ra;(5//3) 



(A.5) 
(A.6) 



where we have replaced r by \S\. 

For points x = Z2 + 5 near the location of the second 
dyon, s = \6\ is small and r = d + 63 + 0{d^/d) is large 
leading to 



and 



47ra,r//3 j sinh(47ru;s//3) 



= e 

^ iwQr/l3 



(A.7) 



e"'~''^ I cosh(47ra;s//3) H — ^ sinh(47ra;s//3) 



(l){x = Z2 + d) 
xix = z2 + 5) 



2d 



Sc0t]l{'i'KLjS/P) + S3 
1 S 

2dsinh{4TrLod/P) 



(A.8) 
(A.9) 
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